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基于 引力 搜索 算法 的 分 数 阶 变 异 时 序 回归 GSA-TSGM(1,1) 模 型 “ 


高 飞 ， 方 海 莲 
(武汉 理工 大 学 理学 院 , AXR 430070) 


摘 要 : 为 了 利用 分 数 阶 累加 算 子 在 灰色 短期 预测 中 的 高 效 性 能 ， 首 次 将 分 数 阶 累加 算 子 引入 变异 时 序 回归 模型 以 期 
取得 更 高 的 预测 精度 。 主要 方法 如 下 : 首先 取 湖 北 省 链子 崖 某 监 测 点 1978 一 1987 年 的 十 年 数据 作为 训练 集 并 使 用 引力 
搜索 算法 确定 最 佳 分 数 阶 累 加 阶 数 ， 而 1988—1993 年 的 六 年 数据 作为 验证 集 验 证 提出 的 模型 ; 其 次 对 比 了 经 典 灰色 模 
型 GM(1,1)、 分 数 阶 累加 灰色 模型 、 变 异 时 序 回 归 模 型 TSGM(1,1) 三 种 灰色 模型 。 结 果 如 下 : 首先 修正 了 陈 西江 等 人 
变异 时 序 回 归 模 型 仿真 时 出 现 的 错误 ， 其 次 表明 了 相 比 于 其 他 的 模型 ， 基 于 引力 搜索 算法 的 分 数 阶 累加 时 序 回归 模型 
在 进行 灰色 长 期 预测 中 具有 较 高 的 预测 精度 。 因 此 ， 通 过 分 数 阶 累 加 算 子 提高 了 灰色 理论 中 长 期 预测 模型 的 精度 ， 为 
灰色 长 期 预测 提供 了 指导 。 

关键 词 : 分 数 阶 累加 算 子 ; 引力 搜索 算法 ; 变异 时 序 回 归 模 型 ; 灰色 长 期 预测 模型 
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Fractional order grey model GSA-TSGM(1,1) with time sequence variation 
regression driven by gravitational search algorithm 


Gao Fei, Fang Hailian 
(School of Sciences, Wuhan University of Technology, Wuhan 430070, China) 


Abstract: In order to use the high performance of fractional accumulated generating operator (FAGO) in the short-term grey 
prediction, in this paper, it firstly added FAGO in the time sequence variation grey model TSGM(1,1) to get higher accuracy. 
The main method organized as follow. Firstly, use the data of 1978 to 1987 from the monitoring station of *Lianziya" mountain 
in Hubei province as training data to optimal FAGO by using gravitational search algorithm and then use the data of 1988-1993 
as verifying data to test the accuracy of the proposed grey model. Secondly, it compared other grey model GM(1,1) , fractional 
accumulated generating GM(1,1) and TSGM(1,1) . The result was that as follow. Firstly, corrected the error in the Chen's article. 
Moreover, it showed that the proposed model in this paper has higher prediction accuracy. Therefore, the novel model improves 
accuracy of the grey theory in long-term prediction by fractional accumulated generating operator and it provides the guide in 
the long-term prediction. 

Key words: fractional accumulated generating operator; gravitational search algorithm; impulse noise; grey model with time 


sequence variation; long-term in grey models 


阶 微分 模型 ， 自 其 提出 后 已 经 广泛 的 解决 了 安全 控制 、 环 境 科 


引言 


Til 


5E. BEUSTETUARAX IRL, G 指 的 是 灰色 ，M 指 的 是 模型 ， 
为 解决 数据 特征 为 不 完整 、 不 明确 的 小 数据 的 预测 问题 ， 第 一 个 “1” 指 的 是 灰 微 分 建 模 方程 的 阶 数 ， 第 二 个 “1” 指 的 


邓 聚 龙 教授 帆 在 1982 年 提出 灰色 模型 理论 ， 其 中 “ 灰 ” 表 示 数 ”是 该 灰色 系统 的 变量 个 数 。 该 模型 有 三 个 特点 : 第 一 ， 可 


据 的 部 分 信息 已 知 , 部 分 信息 未 知 。 与 现 有 的 数据 挖掘 方法 如 ”用 有 限 的 灰 信 息 代 蔡 大 数据 对 事物 未 来 的 发 展 趋势 进行 控 


RK, 分 类 , KKP, 灰色 模型 可 以 利用 有 限 的 、 较 少 的 。 第 三 , 模型 建立 相对 简单 便于 使 用 ; 第 三 , 长 期 预测 效果 不 佳 。 

数据 代 蔡 大 数据 来 实现 预测 。 在 灰色 系统 中 ， 由 于 初始 序列 信 因此， 灰色 理论 作为 一 门 新 兴 的 学 科 并 且 作 为 小 数据 拟 和 

息 的 灰 特 征 ， 常 规 的 拟 合 回归 方法 已 经 难以 对 其 进行 预测 。 经 ”的 有 效 手段 ， 仍 有 值得 学 者 们 进一步 研究 和 改进 的 地 方 。 
灰色 模型 GM(1,1) 是 通过 对 初始 序列 建立 具有 一 个 变量 的 一 随 着 GM(1,1) 模 型 深入 研究 ， 文 献 [4,5] 中 提出 了 不 同 


RI 


以 利 


的 背 
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景 值 会 带 来 不 同 的 系统 误差 。 且 经 


GM (1,1) 模型 中 运 


JE 


小 二 乘法 原理 确定 灰 参 数 本 身 就 会 引入 一 定 的 误差 。 针 对 背景 


值 优 化 ， 谭 等 人 中 首次 构造 
中 经 典 的 背景 值 ， 
增长 序列 也 能 起 到 较 好 的 拟 合 及 预测 结果 。 针 对 
文献 [7,8] 分 别 采 用 扩展 灰 参 数 为 时 间 序 列 和 采用 


出 新 的 背景 值 来 蔡 代 GM(1,1) 模 型 
并 且 新 的 背景 值 不 仅 对 低 增 长 序列 而 且 对 高 
灰 参 数 的 优化 ， 
智能 算法 估计 


灰 参 数 以 减 小 误差 。 而 间接 减少 背景 值 与 灰 参数 带 来 的 误差 也 


是 当时 研究 的 热点 。 文 献 [9,10] 均 是 先 接受 背景 
的 误差 , 但 是 都 通过 残 差 的 原理 进行 分 析 和 优化 。 基 于 文献 [ 


通过 对 时 间 序 列 拟 合 匠 


i 


序 回归 GM(1,1) 的 精度 略 有 提高 。 并 且 
H GM(1,1) 模 型 进行 实例 分 析 时 ， 存 在 一 定 的 


什 与 灰 参 数 带 来 


9], 


陈 等 人 00 认 为 时 序 残 差 GM(1,1) 通 过 原 序列 的 时 序 仍 存 在 一 

的 缺陷 ,于 是 对 1 阶 累 加 生成 算 子 处 理 后 的 序列 时 序 进行 分 析 ， 
3 得 到 变异 时 序 从 而 建 模 提出 变异 时 序 
可 归 GM(1,1) 模 型 。 对比 时 序 残 差 GM(1,1) 的 分 析 结 果 , 变异 
文献 [11] 中 利用 变异 时 
HR ERF 


时 
序 


时 序 回归 模 型 在 预测 未 来 六 年 的 数据 时 可 以 具有 较 好 的 适 ) 


因此 变异 时 序 回归 GM(1,1) 模 型 仍 有 研究 价值 且 其 预测 精 
待 提高 。 
近 些 年 ， 随 着 分 数 阶 微 积分 在 不 同 领域 的 
Wu 等 人 0 首次 将 分 数 阶 微 积分 概念 引入 到 灰 
六 分 数 阶 累 加 灰色 模型 ， 
的 数据 时 , 分数 阶 累加 灰色 模型 可 以 得 到 较 高 的 预测 类 


iu 


泛 应 用 0。 
色 模 型 ， 


we 
HJE o Mao 


等 人 0 将 分 数 阶 累加 算 子 扩展 到 矩阵 形式 ， 
了 方便 。 而 将 平均 相对 误差 作为 目标 函数 3 


F 运 用 


文献 15~17] 均 采 
行 估计 。 

目前 ,灰色 理论 在 研究 人 员 的 不 断 努 力 下 得 以 快速 的 发 
但 是 正 是 由 于 灰色 理论 的 不 断 丰 富 ， 新 问题 也 在 不 断 涌现 。 


为 建 模 的 计算 提供 
DET 
法 进行 估计 这 一 方法 已 成 为 估计 分 数 阶 累 加 阶 数 的 优先 选择 。 
用 粒子 群 优化 算法 对 分 数 阶 累加 算 子 的 阶 数 进 


ES 


文 利用 分 数 阶 累加 算 子 在 预测 方面 的 高 性 能 以 及 目 


模型 进行 改进 ， 将 分 数 阶 累加 算 子 与 变异 时 序 回归 


GMI(1,1 


型 相 结 合 提 出 分 数 阶 累 加 变异 时 序 回归 GM(1,1) 模 型 。 这 样 既 


前 变异 时 序 
习 归 模型 的 预测 精度 有 待 提 高 的 特点 对 变异 时 序 回归 GMCLT) 


) 模 


Hr 
归 GM(1,1) 模 型 的 预测 精度 。 基 于 文献 [18] 中 基准 


累加 算 子 预测 精度 高 的 特点 也 能 提高 变异 时 序 区 
测试 函数 作为 


表明 了 对 于 大 部 分 基准 测试 函 
于 粒子 群 优化 算法 的 。 因 此 ， 本 文采 | 
化 算法 进行 分 数 阶 累 加 算 子 估 计 。 


1 异 时 序 回归 TSGM (1,1) 模型 
1.1 GM(1,1)[3] 

主要 建 模 流程 如 下 : 

GM(1,1) 


输入 :已 有 原 非 负 序 列 X ”= (xod),xz” (2) ,xz (0. 


目标 函数 对 引力 搜索 算法 与 粒子 群 优 化 算法 的 性 能 进行 比较 ， 
数 而 言 ， 引 力 搜索 算法 性 能 是 优 
引力 搜索 智能 启发 式 优 


t. 
RA 
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z< 


a) 规律 化 ， 用 一 阶 累加 算 子 1-AGO 处 理 y” 得 到 X" — 


,4 中 


Gg." Quas e ， 即 a0 - Ys 0 


k z12,..,n 


。 现 "为 单调 递增 序列 。 


b) 建 模 1: 
0.5x (k 2 , k=2,3,..., 


x? (k)+az®(k)=b , KB z” (k) 20.5x (k) + 
n 称 为 背景 值 。 利 | 乘法 估 


最 小 二 
WXS3 ab. HIA B-(-2"Q) e), Y=. O, "OnD)， 


Kp ez0,.D'A JB n-1 的 单位 向 量 。 则 灰 参 数 的 估计 值 


为 P=(a b) =(B'B)'B'Y。 


c) 建 模 2: 将 上 述 估计 的 灰 参数 å, b 代入 灰 微分 建 模 方程 


dx? (t) 
t 


+a Ab 的 时间 o" A 中 


$0( 2 (1) - 5e?» +, r-12,3,.. 


d) 还 原 输出 : 对 得 到 的 序列 名 做 一 阶 累 减 还 原 算 子 


(1-IAGO) BJ $O (r) 2 £P (r1) - $0 (0 , t=2,3,...,1>n 时, 
XQ) 为 预测 值 ， 否 则 为 拟 合 值 。 

以 上 便 是 经 典 灰 色 模 型 的 五 步 建 模 思想 。 
1.2 变异 时 序 回归 GM(1,1) 

可 以 看 出 GM(1,1) 五 步 建 模 中 ， 建 模 的 误差 主要 来 源 于 背 


景 值 和 灰 参 数 的 确立 。 因 此 陈 00 等 人 将 误差 引入 到 时 间 序 列 中 ， 
形成 变异 时 序 。 即 现 先 接受 背景 值 和 灰 参 数 带 来 的 误差 ， 假 设 
根据 GM (1,1) EEA x ^ (n) 没 
主要 建 模 流程 如 下 ; 
变异 时 序 回归 TSGM(1,1) 

a) 时 序 误差 引入 : 假设 况 "w +D 没有 
XO+) » 


误差 。 


ix X B X" (4D) = 


SO =1+6， 


即 g^ C+D= (x9 (1) -byee i E a 5.5 


b) 变异 时 序曲 线 : 上 述 变异 时 序 解 曲 线 为 yo = 
(8960-5 
n - 
à  x"(Dn-2 á 
不 明显 ， 这 里 同样 


tO m 


Qo (r1) - 2) 
x”) -t 
使 用 1 阶 累加 生成 算 子 1-AGO 处 理 o 生成 


o 为 防止 yo 规律 


c) 变异 时 序 回 归 : ”针对 数据 条 目 数 选 取 几 类 常用 的 
2 |n 文献 [12] 中 选取 二 次 多 项 式 f = 
taf ay s 三 次 多 项 式 f =a +at+ at tap UU EE 


9 方程 进行 回归 |， 


S 


B 


f-aytatkagf! af agis SRRIEISTO fat teo A 
Pr212,3..,95 f 519 (2), (3),....,t0 10) 。 数 据 和 曲线 确定 后 
利用 MATLAB 自 带 最 小 二 乘法 原理 多 项 式 曲线 回归 函数 
Popyi(t, 记 可 ， 其 中 心 为 多 项 式 次 数 ， 而 含有 常数 的 宕 函数 为 
非 线 性 曲线 最 小 二 乘法 原理 曲线 拟 合 可 利用 MATLAB 自 带 函 
数 [sgcurvefit 来 实现 。 
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d) 最 佳 变 异 时 序 回归 曲线 : 计算 上 述 备 回归 曲线 的 拟 


ChinaXiv 合 作 期 十 
基于 praster a a DT OX VER (ERIT 


可 见 此 时 7 是 整数 , 根据 分 数 阶 微 积分 组 合 数 的 推广 定义 ， 


合 误差 , 为 防止 过 拟 合 现象 的 发 生 , 文献 [1H] 选 择 拟 合 误差 
为 次 小 误差 的 拟 合 曲 线 作 为 最 佳 变异 时 序 回归 曲线 fo 。 利 用 
fo 曲线 计算 预测 值 ， 最 终 将 得 到 的 序列 做 1 阶 累 减 还 原 (1- 
IAGO) 得 到 zo 的 拟 合 和 预测 值 。 


e) Xe Idi: 代入 条 (ap o (t - 5e 9 +2， 计 算出 


o. IEH 1-IAGO 进行 累 减 还 原 得 到 初始 序列 的 拟 合 及 
预测 值 go 。 


以 上 便 是 文献 [11] 提 出 的 变异 
2 “分数 阶 变 异 时 序 回归 模型 


分 数 阶 累加 生成 算 子 
在 经 典 灰 色 模 型 GM(1,1) 建 模 过 程 中 ， 为 了 使 初始 序列 
x 规律 化 ， 采 用 1 阶 累 加 生成 算 子 对 其 处 理 ， 即 x”) = 


HFE 


We 


模型 的 理 i 


2.1 


2^0 ,采用 箱 阵 的 形式 表示 如 下 : 


i=l 


可 将 + 扩展 成 有 理 数 。 
2.2 ”分数 阶 累 减 生 成 算 子 
同样 地 ， 在 GM(1,1) 建 模 过 程 最 后 需要 对 XO C T ER UR 


还 原 得 到 初始 序列 的 拟 合 预 测 值 。 即 各 "CD 2 590 0 - $00, 
1 
A(O 和 (0) A0 -1 1 ^l A AU 
GODE Qux M= " GODE Qu. (m) (2) 
-1 ls 
Hop ZOBIECADS fü i1. ， 对 角 线 左下 斜 线 位 置 元 素 全 为 


-1 以 及 剩余 位 置 元 素 全 为 0 的 和 矩阵。 显然 ， 根 据 和 矩阵 运算 ， 此 
系数 矩阵 为 4 。 因 此 ，1 阶 累加 生成 算 子 与 1 阶 累 减 还 原 算 子 
互 为 首 算 子 。 同 样 的 ， 根 据 和 矩阵 运算 可 知 ，7 阶 累 减 还 原 系数 
(ny 
累 减 还 原 互 为 送 算 子 。 

2.3 ”分数 阶 变异 时 序 回归 GSA-TSGM(1,1) 模 型 
主要 建 模 流程 如 下 : 


分 数 阶 变异 时 序 回归 GSA-TSGM(1,1) 模 型 
输入 : 已 原 非 负 序列 y” =(x® (D, x” (2) - 3m (n) o 


Id 


矩阵 为 ( 4) = 4 即 x 阶 分 数 阶 累加 生成 算 子 与 7 阶 


a) 规律 化 : 用 


r 阶 累 加 算 子 r-AGO 处 理 


E x" 3 Bl xo 


Z A'XO o 
b) 建 模 1: uq). ar (p) op s RP ug = xk) 
x?(k-1)» z?(k)20.5x"(k)4-0.5x KD k=2,3,...,n。 利 


Ser 


1 
(D (CD Qi 1 1 (0) (0) (0) 
QU DX Dx ADS . (x 0.x Q).sx (7) 
11 1 1d 
0) 
系数 矩阵 设 为 4 ，4 为 下 三 角 元 素 全 为 1 HERE. W r pr 
累加 生成 算 子 为 x (OO = 2370, Bl 
1 
"e 四 1 1 " ; 
Q^ 0,x Qux (e|... Qux Dx S= 
11 1 1 
1 2 
1 1 人 T 
Pus e (x (Dx Qh.sx Q(- 
11 1 1j, 
1 
1 1 [m o o a 
Eu QOD” Oax” M) =A QUPD a" Qua (0) o 
11 1 1 
1 0 0 0 
r 1 0 0 
AKqpa- WP r 1 .0 
: : : . 0 
r(rxl)-G-n-2) r(rxl):-(xn-3) r(rxl):-(r-n-4) 1 
(n-1)! (n-2)! (n-3)! 
; CD ISi 
SA -(a;),.... 则 (a)... » -| j E DIG- )! a 其 证 明 


习 纳 法 参见 文献 [13,14]。 


最 小 二 乘法 估计 灰 参 数 a,p o MS B =h- k) e). 


e=, DV 为 长 度 为 -1 的 单位 


D, A (mm)， 其 中 


灰 参 数 的 估计 值 为 月 = (a P) - army BIY, 


c) 建 模 2: 将 上 述 估计 的 灰 参 数 4 六 代入 灰 微 分 建 模 白化 


线 中 


方程 47^, aunq) cp 的 时 间 响 应 1 $?(Dn- 


(x9 (1) — 2)e gs > {f=1,2,3,...。 


d) 降低 初 值 的 影响 : 对 得 到 的 序列 eo 再 次 做 最 小 二 


Fef 


ib, BH 4 &O (t 1) 2ce?? +d , Bee? e) ， = 


n-1, K'h 


(x^ Q) 0D) 1=1,2,3,..., e 2 (1...,D' 为 长 度 为 


n-] 的 单位 向 量 。 


NX S CR EDS p (è â) = 


(B, By! B, Y, 。 此 时 $O( 1) 2667? cd 


e) 若 此 时 对 eo XETT r 阶 累 减 还 原 ， 则 此 模型 为 文献 [14] 
里 的 分 数 阶 累加 灰色 模型 。 本 文 首 次 对 o 进行 变异 时 序 处 理 
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即 假设 X (rl) 没有 误差 ， 


t” =1+6,， 


) 变异 时 序 


即 太 " +D=XOC+D。 
By s? (16g ads t212,3.... 
线 : 上 述 变异 时 序 解 


线 为 Jo - 


1 T ($9? (1 41) - d) 


a € 


a C 


Br 


HH T , 这 里 同样 使 


(7) 
16 (n 


| 1 阶 累加 生成 算 子 1-AGO 处 理 


d). Aj 


E y 规律 不 


HL hb 
io) 生成 70 。 


HTA 


g 最 佳 


口 


归 曲 线 : 先 使 ) 


对 初始 数据 进行 拟 


FI) 


做 1 fr ERIS JS (1-IAGO) 8 F 
h) 还 原 输出 : 


使 用 
-ATRO 。 
模型 特点 : 


xo 


J2E SEINE E 


归 GM(1,1) 模 型 


合 预 测 得 到 最 佳 时 


曲线 pour o 进行 拟 合 及 预测 得 到 fo 。 
jpo 的 拟 合 和 预测 值 。 


序 回归 上 


最 后 将 得 到 


线 f。 此 时 可 以 
的 序列 


RA g^ aepo ge? 9 qo» WEE go. 


(a) 与 变异 时 序 


态 值 “ 分 数 + ”。 


可 归 模 型 GSA-TSGM(1,1) 
加 算 子 扩展 为 分 数 阶 


I3 


也 就 是 说 变异 时 序 匠 


H 


M— 


日 模型 TSGM(1 


r -IAGO 进行 累 减 还 原 得 到 初始 序列 的 拟 合 及 预测 值 


,1) 相 比 ， 分 数 阶 变异 时 序 


和 持 用 来 处 理 初始 序列 规律 化 的 1 阶 累 
累加 算 子 ， 即 将 原来 的 定 值 “1” 扩充 为 动 


其 实 是 


分 数 阶 变异 时 序 世 


分 数 阶 变异 时 序 回 


n 


归 模 型 GSA-TSG 


Hen 


归 模 型 TSGM(1,1) 的 适用 性 
b) 动 态 值 “分 数 r ”是 


与 初始 序 


分 比 误差 作为 


目标 函数 使 ) 
意味 着 ， 如 果 复 杂 问 题 的 系统 不 同 即 初始 序列 的 训练 集 不 同 ， 
动态 值 “ 分 数 r ”是 不 同 的 ， 即 动态 值 “分数 ”是 


HERRE 


训练 集 决 定 的 
m 


并 控 


型 提高 了 原 有 模型 


的 精度 


c)AE 5E ES] Fr [n 


l 


HA! TSGM(1,1) 


ih MAPE 误差 函数 下 的 。 


归 模 型 TSGM(1,1) 


归 模 型 GSA-TSGM(1,1) 的 特例 形式 。 因 此 ， 


M(1,1) 扩 展 了 原始 变异 时 


Eya El. 


列 中 训练 集 的 平均 绝对 百 
能 算法 优化 决定 的 。 这 就 


初始 序列 
因此 ， 新 模 


是 灰色 长 期 预测 模型 ， 但 


于 模型 精度 


不 


x 


高 所 以 应 用 的 并 不 ) 


在 短期 
ANRI 
灰色 长 


据 测 中 是 提 
[ 算 子 加 入 到 变异 时 序 匠 
期 预测 的 精度 。 


泛 。 而 分 数 阶 累加 算 子 


高 模型 预测 精度 的 
归 模 


强 有 力 工具 。 因 此 ， 将 分 
型 TSGM(1,1) 中 可 以 提高 


3 “引力 搜索 算法 优化 


3.1 目标 函数 

平均 绝对 误差 百分比 (mean absolute percentage error) 定 义 
为 

MAPE = >( xq) - iG) 300) fo -Dxl00% (3) 

文献 [18] 表 明 MAPE 是 最 常用 的 也 是 最 有 效 的 衡量 模型 精 
度 标 准 之 一 。 
3.2 引力 搜索 优化 算法 

和 粒子 群 优 化 算法 相 比 ， 引 力 搜索 算法 也 是 一 种 新 型 启发 
式 智能 优化 算法 。 A Rashedi ADIF 2009 年 提出 了 引力 搜索 
算法 之 后 ， 引 力 搜索 算法 已 经 成 功 地 应 用 于 不 同 领域 的 参数 控 


sio, 


ki 


基于 引力 搜索 算法 的 分 


Chi 
分 数 阶 TR LL GSA-TSGM 
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(1,1) 模型 


在 引力 搜索 算法 里 ， 每 个 对 象 可 看 成 一 个 质点 。 而 每 个 质 
点 在 搜索 空间 里 受 重力 影响 互相 吸引 。 适 应 度 函 数 用 来 代表 每 
个 质点 的 性 能 ， 当 较 轻 的 质点 运动 到 最 佳 位 置 兴 代 过 程 结束 。 
此 时 ， 适 应 度 函 数 已 成 为 全 局 最 佳 适应 度 值 。 
显然 ， 当 引力 搜索 算法 对 分 数 阶 变异 时 序 回 归 模 型 进行 优 
化 时 ， 第 i 个 质点 1 时 刻 适应 度 函 数 如 下 : 
fit.(D) = MAPE(r) (4) 
可 以 看 到 , 适应 度 函 数 是 关于 7 的 一 维 函 数 。 当 7 确定 后 ， 
就 可 以 对 应 一 个 MAPE 值 即 一 个 适应 度 值 。 
引力 搜索 算法 估计 分 数 阶 变异 时 序 回 归 模 型 阶 数 的 流程 图 
如 图 1 所 示 。 利 用 引力 搜索 算法 估计 出 最 佳 阶 数 后 ， 再 代 
入 分 数 阶 变异 时 序 回 归 GM(1,1) 中 求 得 拟 合 和 预测 值 X 09. 
SA xo 
D SENATE. ENT 
DO 并 初始 化 这 N 个 的 位 置 ， 即 在 + 的 上 下 异 范围 内 随机 
化 r 的 值 ; 
b) SERRÉ +=1. 
分 数 阶 变异 时 序 回归 
GM C115 
HRO, 1-12,..N 
即 计算 每 一 次 选 代 时 第 [个 对 象 输出 的 拟 合 预 测 什 
iE fin. PALON 
即 计算 每 一 次 迁 代 时 第 i 个 对 象 也 就 是 r 的 MAPE E 
M 
ERE (HD) Rm xr) 
— 
一 一 ”MAPE 信和 是 否 为 全 局 最 小 值 ? — 
* Yes 
输出 r, MAPEss ) 
图 1 引力 搜索 算法 优化 分 数 阶 变异 时 序 回归 GM(1,1) 流 程 图 
4 ”实例 分 析 
链子 崖 坐落 于 湖北 省 长 江 三峡 沿 岸 ， 以 危 岩 而 成 为 著名 
的 风景 区 。 正 因为 其 地 理 位 置 以 及 危 岸 体 的 存在 ， 安 全 问 


题 成 为 了 众多 地 质 灾 害 专 家 研究 


的 首要 对 象 。 本 文 以 其 中 一 个 


GM(1,1) 的 性 能 。 则 
型 参数 ， 然 后 1988 年 至 1993 年 
索 算法 优化 的 模型 。 并 分 别 给 出 


监测 点 的 观测 数据 作为 实例 来 验证 分 
中 1978-1987 年 的 数据 作为 训练 集 来 得 到 模 


数 阶 灰色 变异 时 序 模型 


的 数据 作为 验证 集 验 证 引力 搜 
GM(1,1) 模 型 ， 分 数 阶 累 加 模 


型 FAGM(1,1) 模 型 ， 变 异 时 序 臣 


归 模 型 GM(1,1) 的 拟 合 预 测 结 


果 ， 与 之 对 比分 析 。 
4.1 修正 文献 [11] 的 错误 

a) 对 初始 序列 用 GM(1,1) 模 
P = (—0.14363,4.8889) 以 及 X? 


接 下 来 首先 修正 文献 [11] 中 出 现 的 错误 。 


型 进行 估计 ， 可 以 得 到 参 数 
= (0.2,5.4883,1 1.5934,18.6414, 


26.7781,36.1716,47.0159,59.5352,73.9881,90.6734)' , 


b) dis 2 jr IRE FF [E 


lH GM(1, 


[o 


1) 模 型 中 的 7 公式 计算 可 得 


= (0,0.8056, 1.6570, 2.5860, 3.8404, 4.9030, 6.0670,6.8924., 
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AEH 


录用 稿 
8.0008,8.9285)' 。 4.2 ”模拟 和 预测 
从 而 1'? = (0,0.8056,2.4626, 5.0486, 8.8890,13.7920, 表 1、2 将 给 出 GM(1,1), 分 数 阶 累加 FAGM(1,1), 变 异 时 序 
19.8590,26.7514,34.7522,43.6807) 。 回归 GM(1,1), 分 数 阶 变 异 时 序 GM(1,1) 的 模拟 值 以 及 预测 值 对 
9 对 1 加 进行 二 次 、 三 次 、 四 次 多 项 式 拟 合 ， 含 有 常数 的 震 比分 析 。 其 中 分 数 阶 累加 FAGM(1,1) 累 加 阶 数 为 >= 0.7060 。 
函数 拟 合 ， 拟 合 曲线 如 下 : 分 数 阶 变异 时 序 GM(1,1) 模 型 的 累加 阶 数 为 
10 20.1085 + 0.10671 + 0.52727 (5) r 21.0365 ,10 的 拟 合 曲线 为 1? —0.4785 — 0.25081 + 0.611360 
1 = 0.4775 — 0.24671 + 0.61111” — 0.00566 (6) —0.00567° , 
1^? = 0.3693 — 0.098471 + 0.55221? + q 仿真 结果 表明 : 
0.003234 刀 一 0.00044131 DDAN RI FAGM(1,1) 模 型 在 进行 长 期 预测 时 ， 精 度 不 
1? 20.570817? +0.172 (8) 如 变异 时 序 回 归 TSGM(1,1) 和 分 数 阶 GSA-TSGM(,1): 
其 中 红色 标注 的 地 方便 是 文献 [11] 中 的 错误 ， 结 果 修 正 会 在 之 b) 变 异 时 序 回 归 TSGM(1,1) 与 分 数 阶 GSA-TSGM(1,1) 的 预 
后 的 表 1 中 给 出 。 测 精度 明显 优 于 GM(1,1) 与 分 数 阶 累加 FAGM(1,1)。 
d) 计 算 以 上 四 种 曲线 拟 合 误差 并 选择 次 小 误差 的 拟 合 曲线 c) Fd 5 可 以 明显 看 到 本 文 提 出 的 GSA-TSGM(1,1) 预 测 误 
作为 最 佳 曲线 拟 合 即 三 次 多 项 式 。 差 是 最 接近 于 0 的 ， 即 预测 精度 最 高 。 
表 1 四 种 灰色 模型 的 拟 合 值 对 比 
实际 值 GM(1,1) 分 数 阶 累 加 FAGM(1,1) 变异 时 序 回 归 TSGM(1,1) 分 数 阶 GSA-TSGM(1,1) 
年 化 
x? /cm XV | em Ro -x? XO | em $o- o| XO | em &o -x XO | em P -x 
1978 0.2 0.2 0 0.2 0.2 0 0.2 0 
1979 4.2 5.2883 1.0883 4.1999 7.6348e-012 4.3702 0.1702 4.3721 0.1721 
1980 5 6.1051 1.1051 5.8055 8.0548e-001 4.1521 0.8479 4.1170 0.8830 
1981 6.2 7.0481 0.8481 7.1535 9.5354e-001 7.7167 1.5167 7.6790 1.4790 
1982 9.8 8.1367 1.6633 8.4540 1.3460 8.8229 0.9771 8.7286 1.0714 
1983 9.8 9.3934 0.4066 9.7788 2.1179e-002 10.0302 0.2302 9.8853 0.0853 
1984 12.6 10.8443 1.7557 11.1657 1.4343 11.3371 1.2629 11.1382 1.4618 
1985 12.5193 2.2193 12.6402 2.3402 12.7398 2.4398 12.4806 2.1806 
1986 15.9 14.4530 1.4470 14.2229 1.6771 14.2320 1.6680 13.9046 1.9954 
1987 15.4 16.6853 1.2853 15.9324 5.3239e-001 15.8042 0.4042 15.4000 1.3838e-010 
MAPE 15.0823% 10.3934% 11.6244% 11.4152% 
表 2 四 种 模型 的 预测 值 对 比 
实际 值 GM(1,1) 分 数 阶 累加 FAGM(1,1) 变异 时 序 回归 TSGM(1,1) 分 数 阶 GSA-TSGM(1,1) 
年 份 
x /cm 2 /cm p" -x® XO | em 2O -40 XO Il em $0 - x0 20 /cm am -x? 
1988 18.1 19.2625 1.1625 17.7868 0.3132 17.4443 0.6557 16.9540 1.1460 
1989 21.3 22.2377 0.9377 19.8046 1.4954 19.1365 2.1635 18.5507 2.7493 
1990 20.1 25.6724 5.5724 22.0047 1.9047 20.8621 0.7621 20.1715 7.1455e-002 
1991 22.0 29.6377 7.6377 24.4077 2.4077 22.5991 0.5991 21.7944 0.2056 
1992 22.6 34.2154 11.6154 27.0355 4.4355 24.3223 1,7223 23.3952 0.7952 
1993 21.4 39.5002 18.1002 29.9117 8.5117 26.0037 4.6037 24.9468 3.5468 
MAPE 34.8734% 14.7619% 8.2382% 6.7702% 
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— Actual values 
16 上 | —*— GM(1,1) 

—$— FAGM( 1) 
14- | —— TSGM( 1) 
—*— GSA-TSGM(1,1) 


DS 


2 四 种 模型 拟 合 


—9— GM(1,1) 

—S— FAGM( 1) 
—PR— 86M(1,1) 
—4— GSA-TSGM(1 1) 


DS 


40 


3 四 种 模型 拟 合 绝对 误差 对 比 


—HI— Actual values 

—*9— GM(1,1) 

asl | 7 FAeMO. O0 
—9— 186M(1,1) 

—$4— GSA-TSGM(1 ,1) 


图 4 四 种 模型 预测 值 


—*— GM(1 1) 

18r  —e— FAGM( ,1) 
—P— TSGM( ,1) 

[| —4€— GSA-TSGM(1 1) 


图 5 四 种 模型 预测 值 绝对 误 


曲线 对 比 


5 ”结束 语 


本 文 主要 解决 了 两 大 问题 


累加 变异 时 序 序列 拟 合 曲线 的 拟 合 


错误 ， 


。 首 先 修 正 了 文献 [11] 中 对 


1 阶 


以 及 变异 时 序 


本 归 


GM(1,1) 对 实例 的 仿真 结果 。 则 
以 及 分 数 阶 累加 灰色 模型 


次 , 本 文 基 


于 变异 时 序 回归 


n 


模型 


首次 将 分 数 阶 累加 算 子 引 入 到 


变异 时 
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分 数 阶 变异 时 序 回 归 GSA-TSG (1,1 


利用 启发 式 智能 算法 进行 估计 。 


归 模 型 中 ， 


序 回归 模型 具有 较 高 的 精度 。 
色 长 期 预测 提供 了 指导 。 
可 以 看 到 ， 在 引力 搜索 算法 对 


阶 数 进行 估计 时 ， 变 异 时 序 回归 的 拟 合 


[分数 阶 


实例 表明 ， 本 文 提 出 的 基于 引力 搜索 算法 的 分 数 阶 变异 时 
同时 为 分 数 阶 累加 算 子 在 进行 灰 


变异 时 序 回归 模型 的 
线 是 在 变异 时 序 回 归 


nu 


模型 进行 仿真 时 提 是 当 


前 确定 了 的 。 但 是 当 
数 是 不 断 变 化 的 ， 最 佳 变异 时 序 回 


归 的 拟 合 


算法 优化 过 程 中 累加 阶 
合 曲 线 可 能 也 会 随 之 


变化 。 因 此 ， 针 对 变异 时 序 回归 的 拟 合 


线 方法 的 选取 仍 值得 


进一步 探讨 。 
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